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I. INTRODUCTION 


It is anticipated that launches associated with the Orbital 
Flight Test (OFT) missions of the Space Shuttle will be conducted 
under wind profile conditions that are less severe than the design 
criteria wind profile. This will require the establishment of tech- 
niques to minimize the probability of exceeding design maximum wind 
loading during ascent. Reduction of wind loading can be achieved by 
wind biasing the ascent trajectory. Ideally, the maximum reduction 
of wind loading would be achieved if the wind profile "seen" by 
the ascending vehicle is known prior to launch. This ideal can 
only be approximated in view of the temporal variability of the 
atmosphere, limitations in available measurement techniques and 
the time lag associated with implementing changes in the ascent 
vehicle wind bias program. However, acceptable wind loading 
can be achieved over most anticipated winds aloft conditions by 





designing a pre-launch wind monitoring plan which provides an 
estimate of in-flight winds within specified error bounds. The 
establishment of the wind monitoring plan will be based, in 
part, on knowledge of the statistics of wind change with respect 


to time. 


This study of wind change over Cape Kennedy, Florida, 
is based on a large sample of winds aloft data (14,610 Rawin- 
sonde profiles) obtained during a fifteen year period. Wind 
change is expressed in terms of component change, unconditional 
and conditional joint distribution of component changes, 
modulus of vector chenge and the joint distribution of wind 


shear component change. 


This report consists of a brief discussion of technical 
background (Section I1]), an analysis of wind change statistics 
(Seetion TIT), a diseussion of conclusions and recommendations 
(Section IV), and listings of the calculated monthly statistics 
of wind change with respect to time at 1 km altitude increments 
from O to 27 km (Appendix). 
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LL. ‘TECHNICAL BACKGROUND 


A. DATA 


Wind change statistics for periods from 12 to 72 hours 
are caleulated from the serially complete (0-27 km) Rawinsonde 
data at To km altitude intervals obtained twice daily during 
the period 1956-70 at Cape Kennedy, Florida. The Rawinsonde 
data obtained four times daily during the period 1962-66 are 
used to verily extension (to time intervals of six hours) of 
theoretical distributions based on the twice daily 1956-70 
data. Sequential Jamsphere wind profile data are used for 


analysis of wind changes for periods less than stx hours. 
B. COORDINATE SYSTEM 


The basde winds aloft data are recorded in terms of 
wind direction, @C and magnitude, W.. The wind vector is 
expressed ino the standard meteorological coordinate system 
in which the dtreetion Crem whieh the wind is blowing is 
measured ino deyrees clockwise from true north. The zonal () 
component, uu, OF the wind veeter is positive for a west (west 
ta Enst > wine Co=P7a" ) wand negative for an east Ceast to west ) 
wind ¢" OO). the meridional component, vy is positive for ia 
south csouth te north) wind (n=180°) and negative for ao norta 
(north to south) wind (a ay, Woand voare obtarned from @ and 
Woaveeord rap te 


u- -Wesin, Os ds 360° (1) 


RS. fe W q’ads, ¢ (2) 


' 


The relation between @& defined above and the 


nogste debined ino the standard mathematical polar form is: 
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ORIGINAL PAGE Is ; 
OF POOR QUALITY .. 
boo time’ The subseript OQ ots used te denete the tnitial value - 
of aovariabte and the subseript To denotes the vartable after 7 
vn elapsed time, At. Thus: e 
‘ : a L 4 ‘ ‘ 
\u uy ty (45 , 
AN. = ee. 5 
\\ uy i (S) ; 
ry : : 
Where. Aucand AY are the components of the wind 3 
change for a specified At. The modulus, R, of the wind bs 
ehanee wrth respeet ta time is given by: e 
R C\udy™ # uAw)e (6) . 
The term wind shear is used exclusively tn this 
report ta deseribe the change of vector wind with respect : 
(oo oa specified vertireal distance below a spectfied altitude. 
The modulus Woo of the vector wind shear ts : 
ri 
‘ ; ey 1% $ 
KR Cu yo e (vi yr (7) ; 
. ; 
] 
Where, a oods the conal wand shear and Vvooois) the 
merrd ream) word sduesar, lt ois eanventional tn discussions 
of wand shear os leulations G0: use the Cerm veetor wtnd shear 
too orepresent the modulus of veetor wind shear, 3 
Yomab and merrvdtenmad wind shear change with respect 
tev Limes are denoted as follows: 
Aw! ue sou! (8) ze 
. t i) . 
AAS = Qe eee (9) : 
1 a) : 
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The modulus of vector wind shear change with respect to 








time is i 
a a eeT| ed 
2 K = Yau) + (Av') (10) 
1 ; 
¥ The means are denoted iy an overbar, the standard 
| deviations and the correlation coifficients are denoted by 
1 o, and R(X,Y), respectively, with X and Y replaced with 
i 
ic the notation appropriate to the variable of interest. 
t D. STATISTICS 
I The wind vector measurements at an initial time and 
f. after an elapsed time are treated in this investigation as 
. a sample from a quadravariate normal distribution defined 
ie by the fourteen statistics listed below: 
MFANS 
a or 7 ee 
Oo’ O° 1’ | | 
f 
é- 
STANDARD DEVIATIONS 
\ 
a OT J 0 
Us) Vy) Uy Vy 
‘ORRFLATION COLFFICIENTS | 
V u u . 
R(UG. 9? , R ( 0’ Dr 
V Vv u v 
R( 0’ 1? i R ( i Se » 
u V u ; 
‘, 
; 
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| 
a 4 : 
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The fourteen statistics of the quadravariate normal 
Q distribution of veetor wind difference with respect to time 
of the five bivariate normal statistics of vector 
nd R Vv 
and R(u,, Vvo)) 


= CONSEST 


: : initial time (u vo. oo ~. AVS 
wind at ano in ( ae he uy Vo 


and the nine statistics involving component differences 


A, Me «J 








= 
which can be calcuiated from the quadravariate statistics 
4 listed above according to the following equations: ' 
MEANS | 
‘ ie ae = Sue a 
| oe 0) C (11) 
| A Mia ee Big = Me NG he? 
a STANDARD DEVIATIONS 
. i 
4 a |. 
a = ; eee A a : , 
Nu Uy, f Uy - Uy, Ny R (uy. Uy) (13) 
bs , 
If 
7% ‘ ° + 3 aS , 47 n . 1 } ’ 
\\ \y Yo ~ Vv, Yq 8 (Vy. Vo) (1-4) an 
Where FON LV) ts the corredtatron coeftieient ; 


Of vatrvabdless N amd ov. 


CORREE VT TON COUEFLCLENTS 


oS Ruy. u,) ae 


Mis gee 
(15) 


VU 





a ns ae rt rele A . 


R (Uy). ‘u) 


Where 


pes arbyt aa tersd from Puauatrenm b.. 





eet eeenetarteetandion sn 
ww & 


mw ote 


ker was 


ae 


nC ee NE i ia ae Se ee ee en | 4 - 5 aw | a | dé wal [ ~ 


5 h 2 


as he 


' 
on —lle 


t « _— we 
Le] ee ania, CIENT aT ctl cS Os. agian 5 


LETT TIFT at 


ae 
ae 





WJ 0 
vy, R (Vo, V1) o 
R (Vo, Av) = ern US 
Tay 
Where oO, is obtained from Equation 14 
o R (vq, u,) _o, RB (uy, v,) 
R (AU, vo) = "I oe OE 
PAu 
J R (u,, ) -9 R (Ug> VQ) 
. Vv 0 Vo 
R (Av, Ug) = ] SRE ae ee eee eee 
ae 
fou, "v4 R (u,> V4) = ar oO. R (u, : Vv? 
+ a a R + 
Uy V1 (uy Vy) "ul °VG R (uy Vo)! 
BE CNS Se 
PAu Cav 
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ITI. ANALYSIS 


INTRODUCTION 


The statistics presented in the appendix of this report 


‘un be useful in the establishment of a basis for certain 


pspects of Space Shuttle Launch planning. A pre-launch wind 


monitoring program may be required to provide data for asséss-~- 
ment or modification of the Space Shuttle wind bias program. 
The development and utilization of the wind monitoring program 
will require knowledge of the magnitude of vector wind change 
with respect to time. The analysis presented in this section 
establishes a theoretical basis for estimation of wind change. 
This is accomplished by comparision of theoretical probability 
distributions, which contain wind change sample statistics 


as parameters (from the appendix of this report), to observed 


probability distributions of wind change. Wind change with 


respect to time is analyzed herein in terms of wind component 
cnange, unconditional and conditional joint distribution of 
wind component change, modulus of vector wind change, and the 


Joint distribution of wind shear component change. 
B. WIND COMPONENT CHANGE WITH RESPECT TO TIME 


The theoretical probatility distribution of wind com- 
ponent change with respect te time ts univariate normal with 
Zero mean and siandard deviation given by Equations 13 and 14; 
the assumption of zero means of component differences is veri- 
fied by the sample statistivs given in the appendix. The 
theoretical normal distribution of component differences can 
bo derived by using either theo standard deviations of component 
ditterences given in the appendix or an estimate which can be 


Obtained Prom the standard deviation of the components if it is 


assumed that: 
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Equations 13 and 14 reduce to 


Oru = a2 or Vi - R(u,, u,) 
ony = ¥2 Wy V1 - R(v,, v4) 


The wind component autocorrelation functions, R(u,, uo) and 





(20) 


R(v,. Vo) can be represented by a negative expontial function 


of time increment, t, i.e., 


R(u,, us) = EXP (-btT) 


R(v,, Vo) = EXP (-cTt) 


where b and c are computed according to 


ut, In R, Cu,» uy) 





b = --~~—-- 


2 


Kxaamples ot the decay of the autocorrelation function 


nt 12 km during January, April and July at Cape Kennedy are 
illustrated in Figure 1; the lines in the figure represent 
the decay rate predicted by Equations 22 and 23. 


Substitution of Equations 22 and 23 into 20 and 21, 


respectively, yields a simple expression for ony and Say 
in terms of o and OL, respectively. 


(22) 


(23) 
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Migure 1. Zonal and Meridional Wind Component Auto- 
correlation at 12 km at Cape Kennedy , 


© Florida (1956-40) 
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On, = V2 9, 1 - EXP (-br) (24) 
Say = 2 0 V1 - EXP (-ctT ) (25) q) 


Equations 24 and 25 indicate that NT and Jay are 
asymptotic to V2 3, andV2 0, for large values of t. Therefore, 
estimates of the extreme value of SAU and O,y are obtained 


by setting 1t equal to ~ in equations 24 and 25. 


The calculated values of b and c for KSC during January, 
April and July are plotted in Figures 2 through 4. The calcu- 
lated and observed values of Tayl™) and OayiT) at 1, 6, 12, 
18 and 24 km during January, April and July are listed in 
Tables 1 through 3. The estimated extreme values of Tau and 
Ony? W2 o, and V2 0, , respectively), are listed at the bottom 
of cach column of calculated values. The comparisons in 
Tables 1 through 3 indicate that Os and nee can be accurately 
estimated by application of Equations 24 and 25, respectively. 
General application of this estimation technique at other 
locations utilizing published statistics of wind component Q 
standard deviations (as in [4] for example) would require a 
more adequate knowledge of the form of the autocorrelation 


lunction than is presently available. 


The cheoretical distribution of wind component differences 
has been derived from sample estimates of Try and Cay and 
Au and Av (;riven in the appendix) for the intervals of 12, 24, 
36 and 48 hours during: January, April and July at 12 km over 
Cape Kennedy; the theoretical normal distributions are plotted 
as straight lines in Figures 5 through 10; the plotted symbols 
represent the observed distributions of Au and Av. It is 
indicat: d that the observed distribution of component changes 
is either accurately or conservately represented by the theoret- 


i, al normal distribution for probabilities from .023 to .977. 
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Table 1. ea eee [Eus. 24, 25) and Observed 
and o during January at Cape Kennedy 
af "1 6, 12, 18 and 24 KM 
"Au Pav 
t(Hours) Calc. Obs. Calc. Obs. 
12 4.09 5.63 4.07 
24 9.63 6.36 4.09 
36 6.72 7.08 4.69 
24 KM 48 7.58 7.56 * 4.67 
60 8.27 8.14 5.06 
72 8.85 8.51 4.99 
oo 12.91 - 
12 Decal 6.60 4.59 4.44 
24 7.06 7.36 5.96 5.40 
36 8.15 8.35 6.75 6.29 
18 KM 48 8.90 8.94 7.24 6.94 
60 9.43 9.49 7.56 7.49 
72 9.82 9.78 7.78 7.75 
oo 11.02 - 8.23 = 
12 7.69 8.70 11.48 9.94 
24 10.51 11.62 14.93 14.20 
36 12.45 13.34 16.92 16.54 
12 KM 48 13.92 14.47 18.17 17.86 
60 15.08 15.18 18.98 18.88 
72 16.02 15.87 19.53 19.36 
o 21.17 - 20.70 = 
12 6.05 6.58 7.53 7.73 
24 8.13 8.45 9.79 10.29 
36 9.49 9.90 11.09 11.48 
6G KM 48 10.47 10.80 11.90 12.18 
60 11.19 11.22 12.44 12.27 
72 11.75 11.76 12.79 12.32 
a 13.86 - 13.55 = 
12 6.92 5.64 7.70 9.82 
24 8.52 8.04 8.63 7.73 
36 9.23 9.31 8.85 8.86 
1 KM 48 9.58 9.70 8.90 9.23 
GO 9.75 9.70 8.92 9.10 
72 9.84 9.62 8.92 8.87 
a 9.93 ~ 8.92 = 
*Validity of Eq. 25 not supported by the data 
at. 24 KM 
14 
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@ Table 2. Culeulated [Bas.2t, 25) and Observed 


rer and oy, During April at Cape Kennedy i, 
Ps 
at i, 6, 2. VR and 24 AM a 
ea 
b| } . 
A AV ri 
( (lours) Cale. Obs . Cale. Obs. 

Le 5 Se | 1.06 2.86 3.12 

ae 3 ea | 2 i hc 3.96 

bo V7 3 a | 3.86 3.90 

241 KM 8 On A ee 1.02 3.905 

va o.8o O.8i titi 4.16 

te Gg.iy i. it 4.15 4.24 

we Tt - 4.20 = 


wre - ee ee es : =~ meee - a ee ae A a oa ee me RS reli MOL 2 he Snel, 


L:! 4.08 bh sn | 4.00 
24 Or O27 9.16 9.01 
36 GO.00 Gol 9.87 9.81 
18 KM 18 Tae TAG CS 3 | 6.34 
uo 7.08 7.80 G.05 G.62 
ie BAT 7.8 OL87 G.78 
11.10 - 7.40 - 
2 LS BH 8.31 a | V8l 
| 1 Sree bia) ae 12.88 
30 1.8 Pk. Fo 1-42.60 Li. 82 
12) KM As | Se | | Ae | | LO | 1.05 
Gu 1 an 16327 re. 9a LST, 
af fg Ve. oe We. 1d tT. 59 7.34 
a Agave - 1.71 - 
}:° ORS yes ty. tnt Yoo 
4 ae: a. ier) re | 
36 er wa ee | 8.57 
ti KM 4S ee 8 O70 SR O18 
aT) LPO LoS POLO Vcct: 38 
1 PP. o LO. 9 eo | 0,33 
ws my. O9 . tO. 31 ~ 
a yo 16 ALYY 18-4 4.69 
. 4 | TS Be) b oO8 OS 
rte 00 S07 u.ag TVS 
1 KM Is SoG Nee 7.02 7.5 
ao SND Soo ies. e000 
oo 8.76 earl 7.380 7.48 
a Ws re be | 





am * 
qa 


al i is a A teat ey va | > ¢ sae ; 
| Lt pL epg ye a oe ae momma , Or i wee : « ~~ ve i it. ‘req 1 ne a ; 





mE 
8 
Table 3. Calculated (Eqs. 24, 25] and Observed - 
OAy and oay during July at Cape Kennedy © 5 
at 1, 6, 12, 18 and 24 KM 7 
| Tau “Ay : 
t (Hours) Calc. Obs. Calc. Obs. ‘ 
12 4.06 4.11 e 
24 3.91 3.59 ‘ 
36 4.18 4.13 : 
24 KM 48 4 4.13 bs 3.57 - 
60 4.53 4.15 t 
72 4.36 3.75 ; 
on ae sf n 
= PSS te ye gS pl ST sy SS ncn GSSOSAEMGN ielPenityren-—s> enast-OaS f 
12 3.17 3.84 
24 2.99 3.14 P 
26 3.50 3.95 ; 
18 KM 48 : 3.66 : 3.63 "s 
60 3.78 4.03 : 
72 3.88 3.79 
e : X 
if 
12 6.76 6.46 5.54 5.49 A 
24 8.97 8.23 7.27 6.90 - 
36 10.34 10.17 8.29 8.27 GQ: 
12 KM 18 11.26 11.08 8.96 8.97 é 
60 11.92 11.98 9.42 9.61 
72 12.39 12.29 9.73 9.85 
w 13.77 10.51 : 
12 4293 3.45 3.68 3.66 
24 4.43 4.21 4.66 4.12 
36 5.13 5.12 5.16 4.94 
KM 48 5.60 5.54 5.45 5.30 
60 5.94 5.95 5.61 5.59 
72 6.19 6.19 5.71 5.69 
~ 6.97 is 5.85 : 
12 3.09 2.95 2.95 3.06 
24 4.09 3.46 3.74 3.37 
36 4.71 4.45 4.14 4,06 
1 KM 418 5.12 4.95 4.36 4.26 
60 5.42 5.51 4.49 4.51 
72 5.62 5.74 4.57 4.49 
0, 6.22 ss 4.68 . 


*Validity of Eqs. 24 and 25 not supported by the 
data oat 18 and 24 KM 
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(1956-70) 


17 





ft 
* # @ 
fo: 

+ 


4 4 
a > 
a 


af ee 


Seer ota 


2, a 
hax 


we Fae, tava 
%, ft  Tak 
+, aE se feat ey hs 


re 


tote gy Re we wav at 
in ett = 
ae ee ind 


- : 
Se # 
‘a 


aif 
reef 
rr 


eam ras ere 
44,8 FR wees Soh 
i eee MEAs 


oO bee 
as a 


Em mae om nemee ineh ire Meet on ER oe ae ar TL RRR, cost be cnn) vie 
7 ’ 
* 
o *T ty : 


7 yt are 
ae 











ye 


‘3 ae r 
aie to * 8 - 24% - 
eae ke sideman sitlerieth Aitidiaaniiny. 6 Oe! gk a itn 


a 
2. hn 


. 
: 

wool ytin ot 
oh 1 OE 








K Gunes} 
eo 


Figure 6. 


it, 





Lae: a. so? q 4 Z 

1y, 7“ * -? of 7 

F 7, é i, Lid ad r { aoe of Kites. | y 1% F te 
.- & es ~ & ae #ar > ors + ur > Mn Be wines ee 


iF 





Ps Ne) ee Nee SEEN eee LN CE LD MME GOR ORME WOR ENY Mme Cem ReR Smee! Vie! | 
eo 6s 1 2 6 10 20 30 46 60 60 70 80 0 069900 (998 08 88 08 8 988 
w[Avcx} GA 9627 


Theoretical (straight lines) and observed 
(plotted points) cumulative probability 
distribution of meridional wind component 
change, Av, with respect to time increment, 
1, during January at 12 km at Cape Kennedy 
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Figure 7. Theoretical (straignt lines) and observed 
(plotted points) cumulative probability 
distribution of zonal wind component 
changes, Au, with respect to time increment, 
1, during April at 12 km at Cape Kennedy 
(1956-70) 
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a Figure 8. Theoretical (straight lines) and observed 
(plotted points) cumulative probability 

rt distribution of meridional wind component 
ss change, Av, with respect to time increment, 


., during April at 12 km at Cape Kennedy 
(1956-70) 
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; Figure 10. Theoretical (straight lines) and observed 
(plotted points) cumulative probability 
distribution of meridional wind component 
change, Av, with respect to time increment, 
t, during July at 12 km at Cape Kennedy 
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Cc. JOINT DISTRIBUTION OF WIND COMPONENT CHANGES WITH 

RESPECT TO TIME 

The joint distribution of zonal and meridional wind 
component changes with respect te time (Au and Av) can be 
upproximated by a bivariate normal distribution. A useful 
property of such a distribution is that an ellipse can be 
calculated which contains the end points of a specified 
percent of vectors having components Au and Av. A detailed 
description of the derivation of probability ellipses and 
plotting, methodology is given by Smith [2). The five para- 
meters of the bivariate normal distribution of Au and Av, 
‘adeulated for each monthly reference period at Cape Kennedy 
at To km altitude intervals from 0 to 27 km are listed in 


the appendix. 


The degree of approximation of the bivariate normal 
distribution to the observed distribution can be evaluated 
by comparison of the observed percentage of vectors which 
are contained within the ellipse to that predicted by the 
ellipse atoa specified probability level. For example, 
for a sample of 1,000 vectors, 950 of the vectors should 
terminate within the 95 percent (theoretical P = .95) ellipse 
ecatealhated from the bivariate statistics of the 1,000 vectors; 
however, a plot of the 1,000 vectors could indicate that only 
AD veetors Cobserved P=.955) terminate within the 95 per- 
cent ellipse. For tllustration on a linear graph comparison 
of the theoretical to the observed P is tiven in terms of 


the parameter as riven by 


=z V- Ln (1-P) (26) 


A comparison of theoretical and observed values of 
a for January, July and April at 12 km for time intervals 
of 12, 24, SG and AR hours is illustrated in Figures 11 thru 13. 
Perfoet agreement between theoretical and observed os is 
represented by a tine drawn from the origin with a slope, B, 
eqauakb to ob. The cateulated Teast squares slopes are given in 
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of Wind Component Changes (Au, Av) with 
Respect to Time at 12 KM During April 
(1956-70) at KSC 
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Figure 13, Observed de as a Function of Theoretical 
\, for a Bivariate Normal Distribution 


of Wind Component Changes (Au, Av) with 
Respect to Time During July (1956-70) 
at 12 KM at KSC 
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The plots indicate an agreement between 
95 (A, < 2.4477). For P >.95 
The interpretation 


the figure legend. 
theory and observation for P < 
the theoretical re exceeds the observed \e° 
of these results is that for extreme probabilities the theoret~ 


ical distributions predict fewer wind change vectors terminating 


outside the ellipse than is observed. These results may have 


to be taken into consideration if engineering application of 
theoretical wind change statistics beyond the 95 percent level 


is required. 

The 95 percent probability ellipses for the joint distri- 
bution of wind component changes with respect to time at 6, 12 
18 and 24 km during January, April and July are illustrated 
in Figure 14; the relatively small changes with respect to 
time during July, the similarities between April and January 
and the large changes at 12 km are clearly illustrated. 
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January, April and July 95 Percent Wind Change 
Ellipses for Time Increments of 12, 24, 36, 48, 60 
and 72 Hours at 6, 123, 18 and 24 KM Over KSC 


28 
ORIGINAL PAGE IS 
OF POOR QUALITY 


a 
we 
» ty 
< 
“? 
ht 
“ 
we 
Py 
4 ~ 
Pe’ 
+a" 
= 
* 
‘J 
* 
<4 
angel 
: 
hws 
“ 
oe 
~ 
* 
eS 
20 | 
a 
ey 
we 


ET 
2 


eae a eer 


re ty Ps : _— " _ ‘ = 
x . a ae HEL ‘ os 2s MRE | A EMRE oy 

J s 3 ot fi ” ee eS] ie eae eet Shae | 

” ts « a7 e 





cw tela 1k BR 


LIT 
‘ + 


he 
Fe 


. - < 
«itr AM ont 


ad . 


Ahh - Meares HE” 


> 








gi i ee 7 om eae ee A 
is he Oa fi es hos ne Be ec 
oa + * a 

; U 

pesce: ie Prt ne See | | 
a on Tria the. nae es ons - 





D. MODULUS OF VECTOR WIND CHANGE WITH RESPECT TO TIME 


If wind changes with respect to time have a distribution 
which is bivariate normal, the modulus R, of the wind change 
vector (defined by Equation 6) has a Rayleigh distribution. 

Since the Rayleigh distribution cannot be integrated in closed 
form, numerical] integration is required to cbtain the cumulative 
probability distribution. Derivation of the Rayleigh distribu- 
Lion, viven the five bivariate normal distribution statistics, 
requires summation involving products of the modified Bessel 
function of the first kind. Smith [2] summarizes the basic 
equations for the Rayleigh distribution derived by Wier [3] and 
oxtended by Yadavalli [4] to include the condition for correlated 
variables. The Rayleigh distribution reduces to the integrable 
classical form if it is assumed that the components of the vector 
wind change are independent and that they have zero means and 
equal standard deviations; the classical Rayleigh probability 


density funetion is 





((R) = 2, EXP (-R?/207) R20 (27) 


Integration of Equation 27 from zero to a specified 
value of R yields the cummulative probability that R < R* where, 


0 (28) 


flv 


Pr {r < Re} = 1 - EXP (-R?/2u’) R 


where oOo = oOo = J 
Au AV 


Since the standard deviation of the component dif- 
ference can be expressed as a function of the standard deviation 
of the components (Equations 24 and 25) it follows that 


2? 
Pr {1 e re} se De i SD Nis ace eee (29) 
40,’ [1-EXP (-kt)) 


| PE nes menee seme ——-— 1 a 





ee den hh ae 


See ee ee ey 





where a. and k correspond to either o, and b or oO. and C 


k 
given in Equations 24 and 25. 

An expression for R given a particular probability, 
Pr (R < R*}, is obtained by solution of Equation 29 to obtain 


R=Vo Ao, Yl - EXP (-kr) (30) 


where A, is derived from Equation 26 denoting Pr [R < R*] 
by P 


The choice of O, = Oy and k =c (from Equation 25) at 
12 km during January, April, and July yields the most accurate 
approximation of the cummulative Rayleigh distribution obtained 
by numerical integration of Equation 28 in Reference 1. A 
comparison of the 99, 95, and 50 percentile modulus of the 
wind change vector with respect to time based on the Rayleigh 
(BRquation 28, Reference 1) and the classical Rayleigh (Equation 
29) is illustrated in Figure 15; the rather good agreement 
indicated for April at 12 km for time intervals from 12 to 72 
hours is attributable to the accuracy of the simplifying 


assumptions described above. 


The remaining question is: How well do these theo- 
retical distributions compare with observed distributions? 
Comparisons of observed and theoretical values of R for time 
intervals of 12, 24, 36 and 48 hours at 12 km during January, 
April and July at KSC are given in Tables 4 through 6; column 
Il of the tables contain R calculated according to the classical 
Rayleigh distribution with o equal to the monthly value of o 
ut 12 km and k equal to the decay constant in the monthly 
exponential least squares fit to the v component autocorrelation 
function (Kquation 23); column I was obtained by numerical 
integration of the Rayleigh distribution. It is indicated 
that the observed cumulative distribution agrees fairly well 
with the theoretical distribution fer probabilities less than 


Vv 


Yo: the observed distribution exceeds the theoretical distri- 
bution for probabilities greater than .95. 
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i. CONDITIONAL VECTOR WIND ELLIPSES 


Prior knowledge that environmental constraints necessary 





lo assure the success of a space vehicle launch will be satis- 


fied implies that there is a capability for prediction of environ- 





montal parameters; the prediction can be based on knowledge of 
condittons prior to Jaunch. With regard to winds aloft, prior 


conditions are typically based on Rawinsonde or Jimsphere wind 
profiles. A typical question that could be posed before launch 
is: Given a measurement of the wind vector 12 hours prior to 
launch at 12 km, will the wind vector at launch time be within 








aac 


(he monthly 95 percent reference month wind ellipse? A question 


of this type can be answered if the distribution of vector 


— «eT 


wind components at an initial time, Ty? and at a future time, 


Ty. can be approximated by a quadravariate normal distribution. 





= Given the components of the vector at Th the conditional 
distribution of the vector wind at T, is bivariate normal. 
Smith [1] describes the derivation of the conditional bi- 
. variate normal distribution and documents the computer program 
© used in this investigation for calculation of these distributions. 
Figures 16 thru 18 illustrate the 95 percent conditional 
bivariate normal distributions at 12 km that have been calcu- 
lated for time increments of 12, 24, 36, 48, 60 and 72 hours 
for the months of January, April and July; five vectors were 
seleeted as given initial conditions for ealculations of the 
conditional ellipses. The components of the vectors are defined 
helow: 
Monthly component means given by Falls [4]. 
2. Maximum zonal wind and the corresponding meridional 
wind from the monthly 95 percent vector wind ellipse. 
3. Minimum zonal wind and the corresponding meridional | 
wind from the monthly 95 percent vector wind ellipse. 
4. Maximum meridional wind and the corresponding zonal | 
wind from the monthly 95 percent vector wind ellipse. 
9. Minimum meridional wind and the corresponding zonal | 
wind from the monthly 95 percent vector wind ellipse. 
: 
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The conditional ellipses illustrated at the center of 
Figures 16 through 18 show that if the observed wind vector has 
components equivalent to the monthly mean components (Condition 
1) then 95 percent of the wind vectors after elapsed times as 
large as 72 hours will fall within the monthly 95 percent ellipse. 
Therefore satisfaction of a launch constraint which states that 
the wind vector must be included within the 95 percent monthly 
Cllipse would be assured for periods as long as 72 hours follow- 
ing an observation of a wind vector having components which cor- 
respond to the monthly means. The conditional ellipses based on 
selection of given wind vectors that terminate on the monthly 95 
percent ollipse (conditions 2 through 5) have a Significant pro- 
portion of their area lying outside the monthly 95 percent el- 
lipse; as the time increment increases this proportion decreases 
but remains significant for a time increment as large as 72 hours. 
This implics that a significant proportion of wind vectors will 
not satisfy a launch constraint based on the 95 percent wind el- 
lipse for periods as long as 72 hours (or longer if these calcu- 
lations are extended) following an observation of a wind vector 


which terminates on the 95 percent ellipse. 


The wind direction characteristics of a wind cllipse 
cain be described in terms of the angles associated with wind 
vectors constructed between the origin and the center of the 


Ctbipse (at the eomponent means) and between the origin and the two 


tangent points to the ellipse. The three vectors constructed 


som: 4 Ae angles v ) ale 
rn this manner and the angles Uys Un. Op. AQ, and Mo are 
tblustrated in Figure 19; the range of wind angles, Vp. is 6 


ta The angles 4 Bie AD and ABs calculated from five 95 


() 
[3° R' 1 
percent conditional ellipse for April at 6, 12, 18 and 24 km 


are btisted tn Table 7. 
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January conditional 95 percent wind ellipses at 
12 km for time inerements of 12, 24, 36, 48, 60 


and 72 hours at Cape Kennedy (1956-70) 


Figure 16. 
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Figure 17. April conditional 95 percent wind ellipses at 
12 km for time increments of 12, 24, 36, 48, 
60 and 72 hours at Cape Kennedy (1956-70) 
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Figure tT8. July conditienal 95 percent wind ellipses at 
12 km for time inerements of 12, 24, 36, 48, 60 
wnd 72 hours at Cape Kennedy (1956-70) 
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F. WIND CHANGES WITH RESPECT TO TIME INCREMENTS LESS 

THAN SIX-HOURS 

The only data suitable for an analysis of wind changes 
aloft at Cape Kennedy for small time increments (<6 hours) 
are the sequential Jimsphere wind profiles obtained during 
the period 1964 thru 1970 [5]. A measurement program which 
began in December 1976 at Cape Kennedy will provide ten 
soundings (six Jimsphere and four Rawinsonde) per day one day 
a week for a 20 week period. These data will be analyzed in 
Phase II of this study. 


Wind changes have been calculated at 6 and 12 km over 
Cape Kennedy from the January, April and July Jimsphere 
sequential runs. The list of dates and number of soundings 
for cach sequential set is given in Table 8. Wind changes 
have been calculated from these data in terms of component 
change (Au, Av) and the modulus, R, of vector change (Eq. 6) 
with respect to time; the calculated Au, Av and R as a function 
of time increment t(denoted by "Delta T'') are illustrated in 
Figures 20 and 21. The wind change data plotted in Figures 
20 and 21 do not line up at exact time intervals because the 
Jimsphere soundings comprising the sequential sets are not 
equally spaced with respect to time. Therefore, calculation 
of wind change statistics utilizing this data set requires the 
use of grouped data. The means and standard deviations of 
component differences for January, April and July at 6 and 12 
km tisted in Table 9 were calculated from data grouped by 1 
hour intervals of t centered at 1 = 1,2...5 hours. The statistics 
do not indicate a strong systematic variation as a function of Tt. 
This is attributed to small sample size and non-uniformity of 
sample size as a function of time increment. Ninety-five 
percent confidence intervals for OAu and Ony? calculated from 
these sample estimates, and theoretical values calculated 
from Equations 24 and 25 are compared in Figures 22 and 23; 
itt is illustrated that in most cases the theoretical values 
are within the 95 percent confidence band. 
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Table So danuary, April and duly 
Sequenbrad Gimsphere Runs 
at Cape Rennedy 
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G. VECTOR WIND SHEAR CHANGE WIEU RESPECT "TO 'PIME 


Vector wind shear change with respect to time can be 
» represented by a bivariate normal distribution, the five 
statistics of the distribution are the means, Au' and Av', 
the standard deviations, Ona! and Digs and the correlation 
coefficient R¢CAu', Avi). Caleulations of these statistics 
for 1 km shear at 12 km supplied by MSFC Space Selences 
Laboratory, were reorganised for uttlivation in this study; 
the statistics for January, April and July are Listed in 
Table 10. The 95 pereent wind shear change ellipses derived 
from these hivariate normal statisties are dllustrated in 
Figure 24. It is indicated that the 995 pereentile shear 
chaage iS largest on January and smallest in July; the 95 
percentile wind shear change is approximately 25 percent 
larger ia January in comparison with April. ‘Phe rather 


close spacing of the ellipses during these months tllustrates 





the fact that wind change is relatively independent of time 
increment for time increments from 12 to 72 hours; therefore, 
most of the 1 km wind shear change at 12 km over a 72 hour 


QO period occurs within the first twelve hours. cr 
The January, April and July 95 pereent to km wind shear 
ellipses at 12 km are also illustrated in Figure 24. It is 


indicated that the 95 pereentite To km wind shear iS smaller 


» a om 


than the 95 percent Tokm wind shear change over time tnerements 


from 12 to 72 hours. 















) Table 10. Bivariate Normal Statistics* of 1 km Vector 
Wind Shear Change with Respect to Time at 
12 km Over Cape Kennedy During January, 
: April and July 
r 1956-67 (Period of Record) f 
Au' Try R(Au',Av' ) Av' Gay’ 
, (m/sec) (m/sec) (m/sec) (m/sec) | 
; = 
. 12 -.01 7.86 1584 02 7.55 
24 .02 8.64 . 2166 . 06 7.84 
i January 36 06 9.15 .2391 .06 7.93 
A8 .03 9.04 . 2364 £13 7.85 
if 60 13 8.76 . 1260 .08 7.67 
ia 1956-70 (Period of Record) 
if 12 .06 5.90 -.0509 =<07 5.62 
| 2A 11 6.31 -.0140 =,10 6.01 
| April 36 .13 6.49 .0459 -.14 5.85 
: 48 .19 6.49 -.0019 ~.15 6.15 
Ee 
y 60 25 6.86 -~.0194 -.18 6.27 
ti 1956-67 (Period of Record) | 
YS : 
if 12 .03 3.89 -.0938 -.02 3.84 
| 24 -,08 4.09 - 0678 00 3.82 . 
ri July 36 12 4,22 -.0385 01 4.06 a 
| AB ~.15 4.14 ~.0405 Ol 4.18 \ 
60 14 4.33 ~.0333 .O1 4.05 | 
FCatlculoted from twice daily Rawinsonde data |: 
is 
{ 
; 
| 
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IV. CONCLUSIONS, REMARKS AND RECOMMENDATIONS 
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The analysis presented in the preceding section for selected O 
months and altitudes illustrates how various theoretical 
distribution functions can be used for calculation of wind change 
with respect to time at Cape Kennedy, Florida. The calculations 
can be made by utilization of the statistics given in the appendix 
for any reference month at 1 km altitude increments from 0 to 27 


km. 





The basic underlying assumption for the calculation of the 
distributions is that the joint distribution of the four variables 
represented by the components of the wind vector at any initial 
time «wand after a specified elapsed time is quadravariate normal. 
If the wind vector is specified at an initial time, then the 
conditional joint distribution of the wind components at a future 
time is bivariate normal. Since each of the variables of the 
quadravariate normal distribution is normal and the difference 
of two normal distributions is normal, it follows that wind 
component change is also normal and the joint distribution of 
zonal and meridional wind change is bivariate normal. The 
modulus of bivariate normally distributed variables has a Rayleigh 
distribution. Therefore, the modulus of vector wind change with 
respect to time is Rayleigh. 


Sample distributions based on reference month Rawinsonde 
data obtained during 1956-70 agree reasonably well with the 
aforementioned theoretical distributions. 


The standard deviation of wind component change with respect 
to time is the only statistic required for determination of the 
theoretical probability distribution (normal with zero mean) of 
wind component change. It has been shown that over a large range 
of altitudes that this statistic can be estimated from wind com~ 
ponent standard deviation and the decay constant of the component 
theoretical autocorrelation function (Figures 2-4). The assumption 
Of exponential decay of the autocorrelation function is reasonably 
accurate in most instances to time increments as large as 60 hours 


a | 7 oe 7 JO ot nner. ee | pave | mom 3 es P Se oe ah ‘ 
: J h to oe 2S daar a | 


Pay 


during January, April and July. The exponential decay model is 
not supported by the autocorrelation data at high altitude during 
January and July (refer to appendix, computation set A, R(X,XP) 
and R(CY,YP)). 


The observed modulus of vector wind change with respect 
to time is systemmatically larger than the predicted modulus 
(Section ITIT1.C.) for probabilities greater than .95. This may 
be attributable to inadequacy of the theory or inaccuracies of 
the data which affect the observed distribution at the extreme 
probabilities. If the theoretical distribution at extreme 
probabilities 1s to be used in engineering applicacions, it 
will be necessary to explain these systemmatic differences. 


Wind change statistics calculated from Jimsphere data for 
smatl time intervals (1 < t <S hours) at Cape Kennedy reveal 
(hat extension of the theoretical calculation of wind component 
standard deviation described above to small time increments is 
valid at 6 and 12 km during January, April and July. A new 
sampling program at Cupe Kennedy which began in December 1976 
will provide six additional Jimsphere runs for each of 20 days 
during one day per week thru April 1977. These data will be 
used in Phase IT of this study in the further analysis of wind 


change for small time intervals. 


SAL iss presently under contract Ccontinuation of NAS8-32226) 


to extend this study of winds aloft temporal variability to include: 


Analysis of vear to year variability 
Establishment of wind change statistics for Vandenberg AFB 
Development and application of a classification ter iwnnique 
for identifiention of homogeneous winds aloft data ets 

e Examination of relations between dynamic stability, 


wind shears and gusts at KSC. 


The final report under the expanded study will be published 


in December 1977, 











- . _——— Se ee 
wt > r 
¢ 7 4 ' a 
“ 
14 


- 





nw 
eng 


! - 7 
ate wep s poe — 
fede te Bluse ee ay ARS a Rab a. eg 


7. 


” « 
- ee soot ales 


«2 


Bak te 


i 

! 

| 
Mi 
{ 

| 
et 
a 





:Or | 


V. REFERENCES 


Smith, O.E.: Vector Wind and Vector Wind Shear Models 
0-27 km Altitude for Cape Kennedy, Florida, and 
Vandenberg AFB, California. NASA TMX-73319, July 1976. 


Weil, H.: The Distribution of Radial Error. Ann. of 
Mathematical Statistics, Vol. 25, 1954, pp. 168-170. 


Yaudavalli, S. V.: On Applications of Some Results 
Related to Bivariate Gaussian Density Distribution 
Functions. Int. J. of Control, ist Series, Vol. 5, 
No. 2, 1967, pp. 191-194. 


Falls, U.W.: Normal Probabilities for Cape Kennedy Wind 
Components - Monthly Reference Periods for all Flight 
Azimuths - Altutudes 0 to 70 kilometers. NASA TMX-64771, 
April 16, 1973. 


Johnson, D. and M. Alexander: Seventy Sequential 
Jimsphere Wind Profile Data Sets for ETR (Cape Kennedy ) 
December 1964 thru July 1970. NASA ES-41, August, 1976. 





Co gte pe eee 8 


“ a 
nn 7, 
? = See 


en 
* aot aie t oe 2 “ 
rd 3 er y 


Pe ee 
af eR. bee coh OL 2 8 ee 


APPENDIX I 


ee pi a 
= € ~~ 7 os. Lae FA x ia > 
ee. Mieka’ 50. seen 


° This appendix contains two sets of reference month quadra- 
variate and conditional bivariate normal statistics of 
variables X, Y, XP and YP, at 1 km intervajJs from 0 to 27 km. 
The statistics were calculated from 15 years (1956-70) of 
twice daily KSC serially complete Rawinsonde data. The 
notution for the variable given in Section II of this report 
differs ‘from the notation established for the computer output 
given herein; the notations are compared in Table I-1. 


TABLE I[-1. NOTATION OF VARIABLES 


COMPUTATION SET 


A B 
Variable Text Computer Text Computer 
(Sect. IT) Out put (Sect.II1) Output 
X Ug u(at T) Ug u(at T) 
y Vo v(at T) Vo v(at T) 
XP Uy u(at T+DT) Uy, ~ Up u(at T+DT) 
= Au -u(at T) 
YP v5 v(at T+DT) Vi - Yo v(at T+DT) 


= Av -v(at T) 


Table I-1 shows that the quadravariate statistics of computa- 
Lion set "A™ are for wind components at an initial time and 
after a specified time increment; the statistics for set "B" 
are for wind components at an initial time and wind component 
change after a specified time increment. The reference month 
quadravariate normal statistics at a particular altitude for 
Six time increments (12, 24, 36, 48, 60 and 72 hours) are 
listed in the Tower left of cach page of computer listing; 
the six sets of conditional bivariate normal statistics cor- 
respondings to the six time inerements are listed in the lower 
right. The data were conditioned on monthly means given 

by Falls [4]. The derivation of the conditional bivariate 
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normal statistics for any other given vector involves re- 
calculation of the conditional means according to equations 
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I-1 and I-2; the standard deviations and correlation co- 
efficients do not have to be recalculated because they are 
independent of the given wind vector. 


| [(Alx,xp) — Rix,yp) Rixp,yp)) (xp* — XB) (,/0,, ) 
5 ixp* =x + “+(Roxyp) ~ Roxxp) Rixp,yp)) (yp* — 9B) (ox/oyp)} 
1~ 2 (1-1) 
[Rixp,yp)] 


[(Rly.xp) ~ Rly,yp) Rixp,yp)) (xp* —%B) (ay/oxp) 


‘+(R(y,yp) — Rly,xp) Rixp,yp)) (yp* — 
Vol yp*=V + evel ny xp Rixp.ypl) (vp" — VB) loy/oyp)] (1-2) 
1 — [R(xp,yp)] 2 


where, X, and Y, are the mean components of the conditional distribution, 
xp* and yp* are the components of the given vector and 


0x, Oy, Syp ANd dyy are equivalent to S.D.x, S.D.y, $.0. xp and 
S.D.yp respectively given in the computer listings. 
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